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ABSTRACT 

A  gaseous  bipropellant  rocket  motor  was  operated  utilizing  ethylene 
and  air.  The  motor  was  of  fixed  geometry  and  had  an  injector  which  pro¬ 
vided  for  controlled  proportions  of  the  total  propellant  flow  to  be  in¬ 
jected  at  various  radial  locations.  Operation  was  accomplished  with  three 
types  of  energy  release  distribution  and  the  resultant  oscillations  ob¬ 
served  and  recorded. 

It  was  determined  that  the  mode  and  severity  of  the  observed  trans¬ 
verse  combustion  pressure  oscillations  were  a  function  of  the  location  of 
the  energy  release  (propellant  injection)  and  the  frequency  of  the  excited 
modes.  From  these  observations,  it  was  concluded  that  a  balance  could  be 
achieved  wherein  two  modes  could  be  excited  with  equal  vigor  and  their 
resultant  interaction  would  produce  stable  operation. 
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INTRODUCTION 

In  the  operation  of  rocket  motors^  variations  in  the  comhustion  pres¬ 
sure  about  the  mean  design  value  may  occur.  These  variations  may  be  random 
or  periodic,  the  random  type  being  termed  rough  burning  and  the  periodic 
being  termed  oscillatory  burning.  Oscillatory  burning  has  been  generally 
divided  into  lov/  frequency,  from  approximately  30  to  500  cycles  per  second, 
and  high  frequency,  from  approximately  1000  to  perhaps  20,000  cycles  per 
second. 

The  low  frequency  oscillations  are  attributed  to  an  interaction  be¬ 
tween  the  dynamic  properties  of  the  combustion  process  and  the  dynamic 
characteristics  of  the  propellant  feed  lines  and  injector  (l)**  The  high 
frequency  oscillations  are  related  to  the  resonant  acoustical  modes  in  the 
combustion  chamber  and  thus  may  be  of  the  longitudinal  mode  or  of  the 
transverse  modes  depending  upon  the  geometry  of  the  chamber  (2). 

The  high  frequency  modes  are  sustained  by  the  energy  released  in  the 
combustion  chamber  during  the  combustion  process.  It  has  been  found  that 
the  location  of  this  energy  release  is  closely  related  to  the  mode  and 
amplitude  of  the  oscillations  that  occur  (5). 

The  object  of  the  series  of  experiments  discussed  herein  is  to  provide 
a  well  defined  energy  release  profile  across  the  injector  face  in  a  rocket 
motor  susceptible  to  treinsverse  oscillations,  in  an  effort  to  determine 

*  Numbers  in  parentheses  are  References  listed  in  th..  Eitli.vgruphy . 
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the  effects  of  these  various  energy  release  profiles  on  the  inherent  in¬ 
stability  of  the  motor. 
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METHOD  OF  INVESTIGATION 

General  Discussion 

An  experimental  investigation  was  conducted  with  the  rocket  motor 
shown  in  Fig.  1.  Premixed  ethylene  and  air  were  the  propellants  utilized. 
The  combustion  chamber  leng'th  was  6  inches  and  the  inside  diameter  was 
7  inches,  the  nozzle  throat  diameter  being  0.5  inches. 

Propellant  flow  was  controlled  during  the  operation  of  the  rocket 
motor  by  regulators  in  the  propellant  feed  lines.  By  adjusting  the  settings 
of  these  regulators,  the  operator  was  able  to  control  the  operating  condi¬ 
tions  of  mean  chamber  pressure  and  equivalence  ratio.* 

Capabilities  of  the  propellant  feed  system  determined  the  ranges  of 
mean  chamber  pressure  and  equivalence  ratio  obtained.  The  range  of  mean 
chamber  pressure  was  from  approximately  30  psia  to  200  psia,  while  that 
of  the  equivalence  ratio  was  from  the  lower  inflammability  limit,  (equiva¬ 
lence  ratio  approximately  equal  to  0-5)  to  approximately  2.0. 

Mean  chamber  pressure,  the  pressure  drops  across  the  propellant 
metering  orifices  and  the  pressures  and  temperatures  upstream  of  these 
orifices  were  automatically  recorded  for  each  datum  point. 

High  frequency  periodic  oscillations  of  the  combustion  pressure  were 
converted  to  an  analogous  voltage  by  a  Photocon-Dynagage  system.  This 
voltage  was  recorded  by  an  Anpex  tape  recorder.  Tne  amplitude  and  frequency 
were  determined  by  comparison  of  the  experimental  signal  with  a  signal  of 

*  Equivalence  ratio  is  defined  as  ohc  ratio  of  the  actual  fuel-oxidizer 
ratio  divided  by  the  stoichiometric  fuel-oxidizer  ratio. 
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FIG.  I  ROCKET  MOTOR  INSTALLATION 
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knovm  amplitude  and  frequency.  More  detailed  information  concerning  the 
apparatus  and  data  reduction  is  given  in  Appendices  C  and  D  and  Reference  5- 

For  each  datum  point,  the  following  information  was  obtained: 

(1)  mean  chamber  pressure, 

(2)  equivalence  ratio, 

(3)  frequency  of  combustion  pressure  oscillations, 

(4)  peak-to-peak  amplitude  of  oscillations,  and 

(5)  mode  of  oscillations. 

For  each  of  the  injection  schemes  described  below,  sufficient  data 
were  obtained  to  determine  an  instability  region  (see  Fig.  4  for  a  typical 
region).  Results  for  each  of  the  injection  schemes  were  then  compared  to 
each  other  and  to  the  applicable  results  of  other  investigators  (5,  1,  ll) • 
Conclusions  were  then  drawn  concerning  the  effects  of  energy  release  loca¬ 
tion,  frequency  and  mode  of  oscillation,  and  the  interdependence  of  these 
effects . 


Method  of  Propellant  Injection 

Three  methods  of  energy  release  were  utilized  in  the  course  of  the 
experimental  work.  These  shall  be  referred  to  as  Energy  Release  Profiles  I, 
II,  and  III,  and  they  are  shown  in  Fig.  2.  Profile  I  has  a  constant  amount 
of  energy  release  as  a  function  of  combustion  chamber  radius,  profile  II 
has  a  decreasing  amount  with  increasing  radius,  while  profile  III  is  such 
that  the  amount  of  energy  released  increases  linearly  with  radius. 

During  the  course  of  this  investigation,  the  assumption  has  been  made 
that  the  energy  release  profile  is  proportional  to  the  amount  of  propellant 
injected  per  unit  of  injector  face  area.  This  assumption  is  justified  by 
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FIG.  2  ENERGY  RELEASE  PROFILES 
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two  considerations.  First,  since  the  propellants  are  premixed  before  in¬ 
jection,  they  burn  (thus  releasing  energy)  very  close  to  the  injector  face, 
the  time  necessary  for  the  physical  mixing  being  longer  than  the  actual 
burning  time  (3).  Second,  since  only  the  transverse  modes  are  under  in¬ 
vestigation,  the  radial  location  of  the  energy  release  is  significant  while 
the  axial  location  is  relatively  unimportant .  This  means  that  the  assumption 
is  valid  unless  the  diffusion  of  the  propellants  from  one  radial  location  to 
another,  admittedly  a  slow  process,  takes  place  more  rapidly  than  the  com¬ 
bustion  process . 

Occurrence  of  only  the  transverse  modes  was  assured  by  selecting  the 
length  of  the  combustion  chamber  below  the  lower  critical  length*  for  the 
ethylene-air  propellant  combination  (2) . 


*  The  lower  critical  length  is  that  length,  for  a  given  propellant  combina¬ 
tion,  below  which  no  oscillations  of  the  longitudinal  mode  will  occur. 

For  ethylene  and  air  it  is  6  inches . 
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EXPERIMENTAL  RESULTS 

Instability  Regions 

Figure  3  depicts  a  typical  instability  region.  An  instability  region 
is  a  map  of  all  the  operating  points  for  a  particular  rocket  motor- 
injector  configuration.  Its  coordinates  are  mean  combustion  chamber  pres¬ 
sure  and  equivalence  ratio.  Since  the  region  contains  all  operating 
pointSj  it  can  be  divided  into  sub-regions  of  stable  and  unstable  opera¬ 
tion.  The  sub-regions  of  unstable  operation  can  be  subdivided  into 
smaller  regions^  each  of  which  comprises  an  area  of  constant  peak-to-peak 
amplitude  of  combustion  pressure  oscillation.  Instability  regions  are  a 
convenient  and  comprehensive  graphical  manner  for  plotting  the  data  ob¬ 
tained  from  the  operation  of  a  rocket  motor  configuration.  The  data  ob¬ 
tained  in  the  experimental  vrork  discussed  herein  will  be  presented  in 
this  form. 


Profile  I 

Figure  3  is  the  instability  region  plotted  from  the  data  obtained 
from  the  experimental  work  employing  profile  I.  Operation  with  profile  I 
resulted  in  a  narrow  region  of  instability  between  equivalence  ratios  of 
0.5  and  0.6.  The  lowest  mean  chamber  pressure  at  which  oscillations  oc¬ 
curred  was  approximately  80  psia. 
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Profile  II 

Operation  with  profile  II  resulted  in  essentially  stable  operation 
over  the  entire  range  of  mean  combustion  chamber  pressure  and  equivalence 
ratio  investigated.  Thus,  no  instability  region  is  presented.  The  range 
of  mean  chamber  pressure  was  from  about  30  psia  to  approximately  200  psia. 
Equivalence  ratios  investigated  ranged  from  the  lower  inflammability  limit 
(equivalence  ratio  approximately  equal  to  O.5)  to  about  2.0. 

In  several  isolated  instances,  low  amplitude  oscillations  vrere  ob¬ 
served;  these  were  of  a  transient  nature  only.  These  observations  were 
random  in  occurrence  and  could  not  be  employed  for  establishing  an  in¬ 
stability  region. 


Profile  III 

Operation  of  the  rocket  motor  system  with  the  injector  having  energy 
release  profile  III  resulted  in  the  instability  regions  presented  in 
Fig.  4.  There  v;as  a  narrow  region  of  instability  at  low  equivalence 
ratios  (from  0.5  to  O.7)  that  was  quite  similar  to  that  of  profile  I. 

For  profile  III  the  lowest  mean  chamber  pressure  at  which  oscillations  oc¬ 
curred  was  approximately  I05  psia,  and  the  observed  oscillations  were 
wealier  than  those  at  similar  operating  points  for  profile  I.  For  example, 
at  an  equivalence  ratio  of  0.55  and  a  mean  chamber  pressure  of  I70  psia, 
profile  III  exhibited  oscillations  of  about  15  psi  while  those  observed 
with  profile  I  were  approximately  10  psi. 

In  addition,  there  was  another  region  located  at  equivalence  ratios 
of  1.2  and  above.  The  lowest  mean  chamber  press'ore  at  which  oscillations 
in  this  region  occurred  was  l40  psia.  At  any  mean  chamber  pressure,  the 


FIG.  4  INSTABILITY  REGIONS  OF  PROFILE 
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oscillations  that  occurred  in  the  lower  region  were  always  of  higher  ampli¬ 
tude  than  those  in  the  upper  region.  For  example,  at  a  mean  chamber  pres¬ 
sure  of  170  psia,  the  maximum  amplitude  of  oscillation  in  the  low  region 
was  15  psi  while  that  in  the  upper  region  was  only  10  psi. 

Modes  of  Oscillation 

Determination  of  Modal  lypes 

The  modes  of  a  transverse  combustion  pressure  oscillation  can  be 
determined  by  either  of  two  methods.  The  first  method  consists  of  stra¬ 
tegically  locating  several  transducers  in  the  combustion  chamber  wall 
(see  Reference  6).  From  the  phase  relationship  of  the  resultant  signals 
obtained  from  those  transducers,  the  mode  of  oscillation  can  be  determined. 

The  second  method  consists  of  computing  a  theoretical  frequency  of 
oscillation  for  the  different  acoustic  modes  and  comparing  this  to  the 
frequency  obtained  in  the  experiments.  The  frequency  of  the  various  modes 
may  be  calculated  according  to  acoustical  theory  (2,  5)* 

In  this  series  of  experiments,  mode  identification  was  accomplished 
by  the  latter  method.  The  following  acoustic  frequencies  are  those  calcu¬ 
lated  for  the  motor  geometry  and  propellant  combination  under  considera¬ 
tion  (2,  5): 
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Table  1 

Theoretical  Acoustical  Freq^uencies 


E.R.  =  1 

E.R.  =  0.5 

First  radial  mode 

6300  cps 

5220  cps 

First  tangential  mode 

3200  cps 

2510  cps 

First  spinning  mode 

3200  cps 

2510  cps 

Second  radial  mode 

11250  cps 

9550  cps 

Second  tangential  mode 

4900  cps 

4l6o  cps 

Second  spinning  mode 

4900  cps 

4l6o  cps 

*  E.R.  is  an  abbreviation  for  equivalence  ratio. 

Profile  I 

For  this  profile,  the  frequency  of  the  oscillations  vras  in  the  range 
from  2400  to  2500  cps.  The  region  of  instability  occurred  at  lovr  equiva¬ 
lence  ratios  (0.5  to  O.t).  Jinc^:  the  theoretical  frequency  for  the  first 
tangential  mode,  in  that  equivalence  ratio  range,  is  approximately  the 
same  value,  2510  cps  (see  Table  l),  it  was  concluded  that  the  observed 
oscillations  vrere  of  the  first  tangential  mode. 

Profile  II 

The  isolated  instances  of  instability  previously  noted  for  this  profile 
were  of  various  modes,  the  first  and  second  tangential,  and  the  first  radial 
all  being  represented.  This  raiidom  occurrence  further  contributed  to  the 
conclusion  that  the  fev;  oscillations  observed  were  transient  in  nature  and 


should  be  disregarded. 


Profile  III 


I'fliile  the  equivalence  ratio  of  the  observed  oscillations  varied  from 
0.5  to  approximately  2.0,  the  frequency  of  the  oscillations  varied  from 
about  2300  cps  at  the  lover  ratios  to  approximately  3OOO  cps  at  the  higher 
values.  By  comparison  of  these  values  with  the  theoretical  frequencies 
given  above  it  was  concluded  that  the  observed  oscillations  of  profile  III 
were  the  first  tangential  mode. 

Peak-to-Peak  Amplitudes 

Rating  of  the  three  profiles  with  respect  to  severity  of  oscillation 
was  accomplished  by  considering  the  peak-to-peak  amplitude  of  the  oscilla¬ 
tions  for  each  of  the  profiles  when  operated  under  identical  conditions  of 
mean  chamber  pressure  and  equivalence  ratio. 

Profile  II  exhibited  only  random,  vreak  oscillations  and  thus  is  con¬ 
sidered  essentially  stable  over  the  entire  operating  range.  In  the  narrow 
region  of  instability  at  low  equivalence  ratios  (O.5  to  O.7),  profile  III 
has  higher  amplitudes  than  profile  I  (I5  psi  as  compared  to  10  psi  at  a 
mean  chamber  pressure  of  17O  psia).  One  exception  is  known,  a  small  portion 
of  the  instability  region  at  chamber  pressures  belov;-  approximately  110  psia. 
Since  profile  III  exhibited  an  unstable  region  at  high  equivalence  ratios 
(above  1.2)  and  also  had  stronger  oscillations  at  common  points  of  in¬ 
stability  as  shown  above,  it  is  considered  more  unstable  than  profile  I. 
Thus,  the  order  of  decreasing  severity  of  unstable  operation  for  the  three 
injection  profiles  would  be  profile  III,  profile  I,  and  profile  II. 


DISCUSSION  OF  RESULTS 


General  Discussion 

Oscillation  Sustaining  Mechanisms 

Considerations  presented  in  Appendix  C  and  elsewhere  (l^  5;  7;  87  ll) 
establish  that  energy  is  required  to  sustain  combustion  pressure  oscilla¬ 
tions.  This  energy  is  dissipated  through  complicated  damping  processes 
(turbulence;  heat  transfer,  and  viscous  effects)  within  the  chamber  as 
well  as  that  energy  lost  by  flow  through  the  nozzle.  The  energy  to  sustain 
these  oscillations  is  obtained  from  the  only  source  available,  the  combus¬ 
tion  process.  Various  mechanisms  (l,  5;  7;  8)  have  been  proposed  by  which 
this  transfer  of  energy  from  the  combustion  process  to  the  oscillating 
gases  is  accomplished.  Most  of  these  are  based  upon  the  considerations 
of  Rayleigh  (4).  Rayleigh  in  I896,  while  discussing  the  maintenance  of 
oscillations  by  heat  addition  to  a  gaseous  medium  observed,  'If  heat  be 
given  to  the  air  at  the  moment  of  greatest  condensation,  or  be  taken  from 
it  at  the  moment  of  greatest  rarefraction,  the  vibration  is  encouraged.'' 

For  any  given  mode  of  oscillation,  there  are  discrete  locations,  called 
pressure  antinodes,  at  which  the  greatest  condensations  and  rarefractions 
occur  (see  Fig.  5)-  For  this  reason,  the  location  of  the  energy  addition 
is  an  influencing  factor  in  determining  which  mode  is  more  likely  to  occur. 

In  the  case  of  the  combustion  chamber  described  herein  the  combustion 
process  adds  heat  (energy)  to  the  gases.  Therefore,  the  interpretation  of 
the  above  is  that  the  oscillations  will  tend  to  be  sustained  if  the  energy 


FIG.  5  MODES  OF  OSCILLATION 
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is  added  at  the  pressure  antinodes  of  the  particular  mode  of  oscillation 
under  consideration. 

The  subject  rocket  combustion  chamber  is  a  right  circular  cylinder 
closed  at  both  ends  (see  Fig.  2).  The  number  of  modes  of  oscillations 
that  might  be  present  are  shown  schematically  in  Figs.  5  and  6.  For 
longitudinal  oscillations,  the  antinodal  locations  would  be  planes  perpen¬ 
dicular  to  the  longitudinal  axis  of  the  chamber  and  in  various  axial  posi¬ 
tions  along  the  chamber  depending  upon  the  harmonic  order.  For  transverse 
oscillations,  the  antinodal  locations  would  be  in  various  radial  positions 
in  the  chamber  (see  Fig.  6). 

Longitudinal  oscillations  were  investigated  by  Davis  (5)  with  injection 
(energy  addition)  in  various  axial  locations,  and  it  was  found  that  injec¬ 
tion  near  the  ends  of  the  chamber  did  in  fact  sustain  the  fundamental 
longitudinal  mode.  Injection  near  the  axial  midpoint,  however,  produced 
relatively  stable  conditions  for  that  mode.  This  is  in  accord  with  the 
Rayleigh  criterion  described  above  (4). 

The  particular  motor  utilized  in  the  series  of  experiments  described 
herein  was  designed  to  be  unstable  only  in  the  transverse  modes  since  a 
6"  chamber  length  was  selected  which  is  below  the  lower  critical  length 
for  the  longitudinal  mode  (2).  To  facilitate  the  examination  and  discus¬ 
sion  of  the  results  obtained,  consider  the  injector  face  as  consisting  of 
three  regions  A,  B,  and  C,  as  sho^vn  in  Fig.  7-  By  comparing  the  regions 
of  Fig.  7  with  the  pressure  distributions  of  Fig.  6,  the  effects  of  the 
energy  released  in  each  of  the  three  regions  can  be  determined  by  applica¬ 
tion  of  Rayleigh's  criterion  (4).  Region  A  corresponds  to  the  location  of 
the  pressure  antinodes  for  oscillations  of  the  first  and  second  radial 


MODE 
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FIG.  6  THEORETICAL  TRANSVERSE  MODES 
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modes.  Tbus^  energy  addition  in  region  A  would  tend  to  produce  oscillations 
of  those  modes.  Both  the  radial  and  the  tangential  modes  have  pressure  anti 
nodes  in  region  C,  Thus,  energy  released  at  that  location  would  he  expected 
to  sustain  radial  and  tangential  oscillations.  Energy  released  in  region  B 
would  not  be  at  the  pressure  antinodal  location  of  either  type  of  oscilla¬ 
tion.  Thus,  it  would  be  expected  to  be  rather  ineffectual  toward  sustain¬ 
ing  them. 


Frequency  Considerations 

Equation  B-3j  derived  in  Appendix  B,  gives  the  qualitative  relation¬ 
ship  between  the  energy  required  to  sustain  an  oscillation  and  the  ampli¬ 
tude  and  frequency  of  that  oscillation. 

E  a  x^  f 
o 

wliere  E  =  energy  required  to  sustain  oscillations, 

x^  =  peak-to-peak  amplitude  of  oscillations,  and 
f  =  frequency  of  oscillation. 

References  (6,  7;  8)  state  that  a  minimum  initial  amplitude  is  neces¬ 
sary  for  producing  sustained  oscillations  and  that  it  is  independent  of 
frequency.  If  then,  energy  is  supplied  to  a  system  (at  the  proper  pres¬ 
sure  antinodal  location)  so  that  two  oscillations  of  unequal  frequency  are 
initiated,  a  lesser  amount  of  energy  will  be  necessary  to  sustain  the  lower 
frequency  oscillation  them  would  be  necessary  to  sustain  the  higher  fre¬ 
quency  mode.  On  the  basis  of  equation  B-3;  it  can  then  be  seen  that  the 
mode  having  the  greater  amplitude  would  be  the  lower  frequency  one.  There¬ 
fore,  it  can  be  concluded  that  the  lovrer  frequency  modes  will  predominate 
when  the  energy  supplied  is  minimal. 


(B-3) 
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For  example,  energy  released  in  region  C  would  tend  to  drive  the 
tangential  modes  as  well  as  the  radial  modes .  The  theoretical  frequency 
values  given  in  Table  1  indicate  that  the  tangential  modes  are  of  lower 
frequency  than  the  radial  modes .  For  this  reason  energy  released  in 
region  C  would  tend  to  excite  and  sustain  tangential  but  not  radial  oscil¬ 
lations.  By  the  same  reasoning,  energy  released  in  any  region  that  would 
tend  to  excite  more  than  one  harmonic  of  the  same  type  of  mode  would  more 
easily  drive  the  lower  order  harmonic. 

Considering  both  location  and  frequency,  the  release  of  energy  in 
region  A  would  drive  the  radial  mode,  the  release  of  energy  in  region  C 
would  drive  the  tangential  mode,  and  that  released  in  region  B  would  be 
rather  ineffectual. 

Discussion  of  Results  of  Profiles  I,  II,  and  III 

General  Discussion 

Several  investigators  in  the  liquid  propellant  rocket  field  (5,  7)  as 
well  as  the  solid  propellant  rocket  field  (11)  have  determined  that,  for  a 
particular  rocket  motor  configuration,  the  location  of  the  energy  release 
greatly  influences  the  amplitude,  mode,  and  frequency  of  the  combustion 
pressure  oscillations.  The  investigation  reported  herein  was  concerned 
with  a  rocket  motor  that  was  capable  of  operation  with  the  three  energy 
release  profiles  presented  in  Fig.  7,  while  the  chamber  geometry  remained 
constant.  It  vras  an  attempt  to  demonstrate  clearly  the  effect  of  injection 
location  on  the  transverse  modes  of  oscillation.  The  following  discussion 
is  divided  into  that  for  profile  I,  profile  II,  and  profile  III. 
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Profile  I 

Figure  3  shows  that  energy  release  profile  I  resulted  in  oscillations 
of  the  tangential  mode.  These  oscillations  occurred  betv/een  the  equivalence 
ratios  of  O.5  and  0.6  and  above  a  mean  chamber  pressure  of  approximately 
80  psia;  their  peak-to-peaik  amplitude  increasing  at  higher  mean  chamber 
pressures  and  getting  as  high  as  15  psi  at  a  chamber  pressure  of  I90  psia. 

Profile  I  has  an  equal  amount  of  energj'-  released  per  unit  area  in  each 
of  the  three  regions;  A,  B,  and  C  (see  Fig.  7)-  Thus^  by  the  location 
considerations  given  above,  this  injection  pattern  vrill  supply  equal 
amounts  of  driving  energy  to  the  tangential  and  the  radial  modes.  As  a 
result,  the  relative  magnitude  of  the  frequencies  of  the  two  simultaneously 
driven  modes  become  significant.  Since  the  tangential  mode  is  of  lower 
frequency,  the  energy  supplied  to  it  is  more  efficient  for  producing  oscil¬ 
lations  than  that  made  available  to  the  radial  mode.  Thus,  the  tangential 
mode  is  predominant. 


Profile  II 

Operation  with  profile  II  resulted  in  essentially  stable  operation 
over  the  entire  range  of  mean  chamber  pressures  and  equivalence  ratios  in¬ 
vestigated  even  though  energy  is  released  in  locations  v/ell  suited  for 
sustaining  transverse  node  oscillaticns .  That  is,  region  A  receives  the 
largest  amount  of  the  total  energy  released,  and  this  energy  acts  to  drive 
the  radial  mode  of  oscillation.  The  somewhat  smaller  amount  of  energy- 
released  in  region  C  tends  to  sustain  the  tangential  mode.  The  ability  of 
either  of  the  two  modes  to  become  predominant  is  dependent  upon  the  fre¬ 
quency  of  the  mode  as  vrell  as  the  amount  of  energy  that  it  receives. 
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Profile  II  balances  these  effects,  providing  a  relatively  large  amount  of 
energy  to  the  higher  frequency  radial  mode. 

Thus,  the  distribution  of  energy  as  described  above  causes  what  can  be 
termed  an  intermodal  conflict  between  the  tangential  mode  and  the  radial 
mode.  That  is,  the  pressure  and  velocity  variations  which  occur  for  these 
modes  are  not  compatible  spatially  (see  Fig.  6).  Moreover,  since  the  fre¬ 
quencies  of  these  two  modes  are  unequal  (see  Table  l),  these  pressure  and 
velocity  variations  are  not  in  phase  and  interfere  with  each  other.  Thus, 
this  is  termed  an  intermodal  conflict  which  for  this  particular  injection 
system  resulted  in  stable  operation.  Obtaining  such  an  intermodal  con¬ 
flict  may  become  an  important  means  whereby  liquid  propellant  rocket  motors 
can  be  stabilized. 


Profile  III 

Figure  4  shows  that  operation  the  rocket  motor  using  profile  III 
resulted  in  a  narrow  region  of  unstable  operation  centered  at  an  equiva¬ 
lence  ratio  of  0.6  as  vrell  as  another  region  at  higher  equivalence  ratios 
(above  1.2).  The  mean  chamber  pressures  at  which  unstable  operation  oc¬ 
curred  was  above  approximately  l40  psia.  The  observed  oscillations  were 
of  the  tangential  mode. 

These  results  were  obtained  with  an  energy  release  profile  that  pro¬ 
vided  a  large  amount  of  energy  release  in  region  C  and  a  somewhat  smaller 
amount  in  region  A.  This  situation  is  the  reverse  of  that  presented  by 
profile  II  which  achieved  a  balance  between  the  effects  of  location  and 
frequency.  Profile  III,  however,  supplies  the  lower  frequency  mode  (tan¬ 
gential)  with  the  majority  of  the  energy  available,  and  therefore,  this 
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arrangement  is  more  favorable  for  driving  tangential  oscillations  than  was 
profile  I  or  profile  II. 

f 

Comparison  of  Profiles  I,  11^  and  III 
Comparison  of  the  results  of  profile  III  with  profile  I  disclosed 
that  the  tangential  mode  is  driven  more  strongly  by  profile  III  than  pro¬ 
file  I,  while  the  radial  mode  is  driven  less.  This  comparison  would  indi¬ 
cate  that  the  predominant  tangential  mode  of  oscillationj  which  is  common 
to  both;  should  be  stronger  for  profile  III. 

That  this  is  the  case  can  be  determined  by  comparison  of  Figs.  3  and 
4.  In  general  the  oscillations  for  a  given  operating  point  for  profile  III 
were  stronger  than  those  for  similar  conditions  with  profile  I.  Also,  os¬ 
cillations  observed  with  profile  III  occurred  over  a  larger  portion  of  the 
operating  range  than  did  those  of  profile  I. 

Thus  the  experimental  results  verify  the  predictions  based  on  con¬ 
sideration  of  the  energy  release  location  and  the  frequencies  of  the 
simultaneously  excited  modes. 

Comparison  of  Results  with  Those  of  Other  Investigators 
Davis  (5)  conducted  i  .series  cf  experiments  vdth  the  same  parameters 
of  chamber  geometry  and  propelleint  combination  maintained  constant.  The 
injectors  employed  in  that  program  are  shown  schematically  in  Fig.  7  as 
injectors  IV,  V,  and  VI.  Each  of  these  three  had  one  ring  of  injection 
holes,  located  at  diameters  of  2",  4",  and  6"  respectively,  thereby  cor¬ 
responding  to  injection  purely  in  regions  A,  B,  and  C. 

The  results  obtained  by  (5)  are  in  exact  agreement  with  the  observa¬ 
tion  of  (4).  Since  these  injectors  did  not  excite  two  modes  simultaneously, 
the  frequency  considerations  previously  discussed  were  not  involved. 
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Injector  IV,  injection  on  the  2"  diameter  circle,  produced  instability 
in  the  radial  mode.  Injector  V,  injection  on  the  4"  diameter  diameter 
circle  resulted  in  stable  operation  over  the  entire  range  of  mean  chamber 
pressure  and  equivalence  ratio  investigated.  Injector  VI,  injection  on 
the  6"  diameter  circle  resulted  in  tangential  oscillations. 

Qualitative  comparison  between  the  results  of  the  6  injection  profiles 
presented  in  Fig.  7  is  possible  because  the  parameters  given  in  Fig.  1 
were  maintained  constant  for  all  six  profiles . 

To  illustrate  the  relative  instability,  the  maximum  amplitudes  of  os¬ 
cillation  observed  at  a  mean  chamber  pressure  of  I50  psia  are  given  below: 

Profile  Amplitude  Mode 

(psi) 


I 

II 

III 

IV 

V 

VI 


7 

1st  tangential 

0 

11 

1st  tangential 

13* 

1st  r'idial* 

0* 

40* 

1st  tangential* 

*  From  Reference  5* 


Profiles  IV  and  VI,  which  each  contribute  energy  to  only  one  mode, 
produced  the  strongest  oscillations.  Profile  VI  drove  a  lower  frequency 
mode  than  profile  IV  and  thus  produced  stronger  oscillations  than  did 
profile  IV. 

Profiles  I  and  III  shared  their  available  energy  between  two  modes 


and  thus  produced  weaker  oscillations  than  profiles  IV  and  VI.  Profile  III 
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exhibited  stronger  oscillations  than  profile  I  because  its  energy  release 
pattern  more  heavily  favored  the  predominant  tangential  mode. 

Observations  indicate  that  profiles  II  and  V  produced  stable  operation. 
The  reasons  for  this  behavior  were  not  the  same  in  both  cases.  Profile  V 
was  stable  because  its  energy  was  released  in  a  location  that  did  not 
readily  excite  either  tangential  or  radial  modes  (pressure  node,  see 
Ref.  4).  Profile  II,  however,  balanced  the  effects  of  location  and  fre¬ 
quency  to  produce  stable  operation. 

Pickford  and  Peoples  (7)  discuss  the  fact  that  the  combustion  process 
provides  the  energy  which  drives  and  sustains  combustion  pressure  oscilla¬ 
tions.  They  indicate  that  the  phase  relationship  between  the  energy  ad¬ 
dition  and  the  action  of  the  pressure  wave  itself  is  extremely  important. 

The  phase  relation  between  the  energy  addition  is  dependent,  in  part,  on 
the  location  of  the  energy  release  within  the  combustion  chamber.  This 
consideration  is  acknowledged  by  (7)  when  they  conclude,  while  discussing 
experimental  observations,  that  differences  in  propellant  spatial  distribu¬ 
tion  can  cause  variations  in  the  inherent  stability  of  a  rocket  motor. 

The  observations  and  conclusions  presented  herein  are  in  agreement 
with  the  work  of  (j).  They  are  an  attempt  to  more  clearly  establish  the 
nature  of  the  differences  in  spatial  distribution  mentioned  by  (j).  It 
has  been  observed,  however,  that  the  detrimental  effects  of  energy  release 
location  can  be  countered  by  the  intermodal  action  observed  during  opera¬ 
tion  with  profile  II. 

Experimental  v/ork  concerned  with  the  location  of  energy  release  using 
solid  propellants  has  been  accomplished  by  Price  (ll).  His  program  con¬ 
sisted  of  constructing  solid  propellant  motors  with  charges  made  partly  of 
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dummy  propellant.  The  location  of  the  turning  surface  was  moved  gradually 
closer  to  the  pressure  antinode  of  the  unstable  acoustic  mode.  The  ampli¬ 
tude  of  the  observed  oscillations  was  determined  to  be  a  function  of  the 
location  of  the  burning  surface  (energy  release),  the  amplitude  being  the 
greatest  when  this  location  was  coincident  with  the  pressure  antinodes. 

The  results  of  (..l)  are  in  exact  agreement  with  those  of  (5)*  Both 
investigators  were  concerned  with  the  location  of  the  energy  release  with 
respect  to  the  pressure  antinode  of  a  singular  mode  of  oscillation.  Re¬ 
sults  of  the  subject  experimental  investigation  agree  with  those  of  (5) 
and  (11.)  with  regard  to  excitation  of  a  singular  mode.  They  also  include 
consideration  of  the  frequencies  of  oscillation  if  more  modes  than  cr.e  are 
excited  simultaneously. 

The  observations  of  (5),  (?)>  and  (ll),  and  the  results  presented 
herein  are  ail  in  agreement  with  the  considerations  of  (i+).  While  these 
results  are  significant,  the  other  variables  such  as  geometry,  fuel  com¬ 
position,  mixture  ratio  and  chamber  pressure  should  not  be  neglected. 


CONCLUSIONS 


From  the  results  of  the  experiments  described  herein  as  well  as  the 
results  and  observations  of  other  investigators  (5,  7;  H)  i't  may  be  con¬ 
cluded  that  the  mode  and  amplitude  of  transverse  combustion  pressure  oscil¬ 
lations  are  very  closely  related  to  the  radial  energy  release  variations  of 
the  injection  system. 

Additional  conclusions  were  (l)  Injection  at  the  pressure  antinodal 
location  of  a  possible  mode  of  oscillation  will  tend  to  excite  that  mode; 

(2)  If  injection  under  similar  operating  conditions  in  two  similar  rocket 
motors  is  accomplished  in  a  manner  which  excites  one  mode  in  the  first 
and  a  different  mode  in  the  second;  the  lower  frequency  mode  will  then  ex¬ 
hibit  the  larger  peak-to-peak  amplitude  oscillations,  and  (3)  If  more  than 
one  mode  of  oscillation  is  excited  simultaneously,  the  mode  and  amplitude 
of  the  predominant  oscillation  will  be  determined  by  the  frequency  of  each 
mode  as  well  as  the  amount  of  energy  supplied  to  each  mode. 

Stable  operation  may  be  achieved  by  the  following  methods:  (l)  The 
propellants  can  be  injected  at  radial  locations  that  do  not  coincide  with 
the  pressure  antinodal  locations  of  possible  modes  of  oscillation  and 
(2)  The  propellants  can  be  injected  in  a  manner  that  will  drive  two  modes 
of  approximately  equal  strength  and  the  interaction  between  them  will 


produce  stable  combustion. 
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APPENDIX  A 

NOTATION 

Explanation 

area 

viscous  damping  coefficient 
compressibility  factor 
correction  factor 
cycles  per  second 
energy 

equivalence  ratio 
frequency 
flow  factor 
force 

number  of  holes 
spring  constant 
orifice  constant 
mass 

mean  static  pressure 

mean  chamber  pressure 

pounds  per  square  inch^  absolute 

degrees  Rankine 

time 

temperature,  absolute 


Units 

•  V.2 

inch 

Ib^-sec  per  inch 
none 

none 

cycles  per  sec 

ft-lb^ 

none 

cycles  per  sec 
none 

none 

Ib^per  inch 

none 

slug 

2 

Ib^per  inch 

2 

Ib^per  inch 

2 

Ib^per  inch 
degrees  Rankine 
sec 


degrees  Rankine 


Symbol 

Explanation 

Units 

w 

flow  rate 

lb  per  sec 

X 

displacement 

inch 

X 

dx/dt 

inch  per  sec 

•• 

X 

d^x/dt^ 

inch  per  sec' 

X 

0 

peak  amplitude  of  oscillation 

inch 

Y 

expansion  factor 

none 

t 

specific  weight 

lb  per  foot 
m 

n 

3.14159 

none 

(JU 

angular  frequency 

rad  ner  sec 

3^ 


APPENDIX  B 

FUNCTIONAL  RELATIONSHIP  FOR  SUSTAINING  ENERGY 

In  the  discussion  of  the  combustion  pressure  oscillations  observed 
in  this  experimental  investigation,  the  effect  of  frequency  upon  the 
energy  required  to  sustain  the  oscillations  is  introduced.  There  is  no 
relationship  available  in  the  literature  that  relates  these  two  variables. 
This  situation  is  due  to  the  nonlinear  nature  of  the  damping  processes 
and  the  fact  that  the  nature  of  the  damping  processes  (turbulence,  radia¬ 
tion,  viscosity)  is  not  clearly  understood. 

Consider  the  situation  present  in  a  rocket  motor  operating  with  high 
frequency  combustion  pressure  oscillations.  Within  the  chamber  there 
exists  an  oscillating  system.  That  is,  the  gas  particles  in  the  combustion 
chamber  have  an  oscillatory  motion  superimposed  upon  the  steady  flow  mo¬ 
tion.  The  unbalanced  pressure  within  the  chamber  is  the  driving  force 
which  causes  the  particles  to  oscillate.  This  driving  force  acts  against 
the  inertial  forces  which  are  due  to  the  momentum  of  the  particles.  More¬ 
over,  the  oscillations  will  die  out  if  energy  is  not  supplied  to  them; 
thus,  some  sort  of  damping  forces  are  present. 

Basically  then,  three  elements  exist  which  are  analogous  to  the  three 
elements  of  a  simple  lumped  parameter,  mass,  spring,  and  damper  system. 

The  gases  themselves  constitute  the  mass,  the  unbalanced  pressure  is 
analogous  to  the  spring  and  the  damping  processes  mentioned  above  are 
similar  to  those  processes  which  occur  in  the  viscous  damper.  Certainly 


the  analogy  is  not  very  strong  as  the  character  of  each  of  the  elements 
vary  between  the  two  systems.  For  this  reason  caution  must  be  used  when 
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drawing  conclusions  based  upon  it.  However,  the  analogy  will  be  used  to 
help  establish  a  functional  relationship  which  successfully  aids  in  the 
explanation  of  combustion  pressure  oscillations. 

The  equation  of  motion  of  the  simple  lumped  parameter  mechanical 
system  is  (lO): 

mk’  +  cx  +  kx  =  F  cos  uJ  t  (B-1) 

0 

where  m  =  mass, 

X  =  second  time  derivative  of  x, 

c  =  viscous  damping  coefficient, 

X  =  first  time  derivative  of  x, 

k  =  spring  constant, 

x  =  displacement,  and 

F^cosajt  =  sinusoidal  driving  force  applied  to  m. 

The  amount  of  energy  put  into  such  a  system  may  be  obtained  by  taking 
the  product  of  the  force  applied  and  the  displacement  of  the  force  and 
then  integrating  the  product  over  a  time  interval  equivalent  to  one  cycle. 
The  resulting  equation  is  (lO) : 

E  =  -  2jt^  c  x^  f  (B-2) 

0 

where  E  =  energy  transferred, 

c  =  viscous  damping  coefficient, 
x^  =  peak  displacement, 
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f  =  frequency  of  oscillation, 

and  where  the  minus  sign  merely  indicates  that  the  energy  is  being  put  into 
the  system,  obviously  an  equal  amount  is  dissipated. 

Due  to  the  weakness  of  the  analogy,  equation  B-2  cannot  be  taken  as 
the  equation  which  describes  the  energy  addition  to  an  oscillating  gaseous 
medium,  but  on  the  basis  of  the  fact  that  both  systems  are  forced  oscil¬ 
lating  systems,  the  following  functional  relationship  is  applicable  to 
either: 

E  a  f  (B-3) 

where  the  terms  are  the  same  as  they  were  in  equation  B-2.  Because  equa¬ 
tion  B-3  is  a  functional  relationship,  the  term  represents  either  the 
peak  displacement  from  the  midposition  or  the  peaJ:-to-peak  amplitude, 
these  two  displacements  differing  by  a  factor  of  2. 

Equation  B-3  describes  the  manner  in  which  the  energy  necessary  to 
sustain  oscillations  in  a  rocket  combustion  chamber  is  dependent  upon  the 
amplitude  and  frequency  of  these  oscillations.  It  shows  that  the  energy 
required  to  sustain  a  larger  amplitude  of  oscillation  varies  with  the 
square  of  the  peak-to-peak  amplitude  while  this  energy  requirement  is  only 
a  first  power  function  of  frequency. 

When  considering  the  oscillations  in  a  rocket  combustion  chamber, 
several  different  modes  of  unequal  but  rather  well  defined  frequencies  are 
possible.  An  interesting  situation  arises  when  the  amount  of  energy  made 
available  to  drive  an  existing  oscillation  is  increased.  Two  possibilities 
exist,  the  amplitude  may  increase  or  the  frequency  may  shift  to  that  of  a 
different  mode  or  perhaps  a  higher  harmonic  of  the  existing  mode.  In  the 
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event  of  a  mode  shift  the  amplitude  may  change  slightly  to  satisfy  the 
function  relationship  of  equation  B-3- 
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APPENDIX  C 

DESCRIPTION  OF  APPARATUS 

General  Description 

The  test  cell^  control  rooni,  propellant  feed  system  and  instrumentation 
are  fully  described  in  Reference  (5)- 

Figure  8  is  a  cross  sectional  view  of  the  rocket  motor  utilized  in  this 
series  of  experiments.  The  motor  consists  mainly  of  4  components;  the  pre¬ 
mixer,  the  injector,  the  combustion  charaber,  and  the  nozzle. 

The  Combustion  Chamber 

The  combustion  chamber  is  an  uncooled  right  circular  cylinder  with  a 
7"  inside  diameter  and  6"  axial  length.  It  is  constructed  of  mild  steel 
with  mounting  flanges  welded  to  each  end  to  receive  the  injector  plate 
and  the  nozzle  plate.  Tapped  holes  are  provided  in  the  chamber  wall  for 
the  transducers,  mean  chamber  pressure  tap  and  post  run  chamber  coolant. 

The  Nozzle  Plate 

The  nozzle  plate  was  a  mild  steel  r]ate  v.'lth  a  DeLaval  nozzle  with  a 
0.5"  throat  diameter  v.-elded  to  it.  Ibe  converging-diverging  portion  of 
the  nozzle  was  constructed  of  stainl  ss  steel.  The  converging  and  throat 
portion  tos  convecti^clv  ccolel  with  ’./ater  u.;id  the  diverging  portion  was 
cooled  by  e.xternal  cooling  tvibes.  A.i  a'lto-.otive  typo  spark  plug  was  in¬ 
stalled  in  the  nozzle  plate  to  provide  igr.iticn. 
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The  Premixer 

The  premixer  was  a  stainless  steel  chamber  filled  with  two  hundred 
and  twelve  l/2"  diameter  brass  balls  as  shown  in  Fig.  8.  The  fuel  and 
oxidizer  enter  and  pass  through  the  packed  balls ^  thus  being  thoroughly 
mixed.  The  mixed  propellants  then  flow  through  seven  lines  to  the  seven 
concentric  rings  of  injection  holes  in  the  injector  plate. 

The  Injector  Plate 

The  injector  plate,  shown  in  Fig.  9,  was  convectively  cooled  on  its 
exterior  surface.  The  injection  hole  pattern  shown  in  Fig.  10  consists  of 
seven  concentric  rings  of  holes,  spaced  such  that  each  ring  was  centrally 
located  in  one  of  the  circular  strips  of  injector  face  shown  in  Fig.  10. 
These  rings  are  denoted  as  ring  numbers  1,  2,  3;  5;  8,  and  7-  Kie  areas 

are  in  the  ratios  of  1,  3>  5;  Ij  9,  H;  and  13  and  thus  the  rings  have 
3)  9)  1)  9,  and  13  holes  respectively.  The  holes  were  arranged 

radially  at  random  in  an  effort  to  avoid  local  concentrations  of  holes  and 
thus  energy  release . 

Each  ring  of  injection  holes  issues  from  a  manifold  machined  into  the 
upstream  side  of  the  injector  plate.  Bosses  welded  to  these  manifolds 
receive  the  lines  from  the  premixer.  The  two  outermost  manifolds  each 
have  two  bosses  to  reduce  the  possibility  of  a  sizeable  pressure  drop 
within  these  manifolds  causing  an  uneven  distribution  of  propellant  to 
each  hole. 


Flow  Metering  Orifices 

The  premixer  and  the  injector  were  connected  by  seven  l/U"  stainless 
steel  lines.  Each  line  contained  a  small  brass  orifice  (shown  in  Fig.  8) 


FIG.  9  INJECTOR  PLATE 


INJECTION  HOLE  PATTERN 
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at  the  premixer  end,,  that  became  choked  during  operation.  Since  the 
orifices  were  choked,  the  fraction  of  the  total  propellant  flow  that 
passed  through  each  one  was  a  function  of  the  upstream  conditions,  which 
were  common  to  all  orifices.  Thus,  the  flow  rate  became  a  function  of 
the  size  of  the  orifice  only. 

In  order  to  eliminate  the  possibility  of  various  discharge  coeffi¬ 
cients  affecting  the  selection  of  orifices  on  an  area  basis,  a  large  number 
of  orifices  were  constructed  and  checked  with  the  air  flow  rate  necessary 
to  cause  choking  while  maintaining  constant  upstream  conditions .  The 
values  of  flow  rate  thus  acquired  were  then  used  to  determine  which  com¬ 
bination  of  orifices  would  most  closely  develop  the  desired  profiles. 

For  each  profile,  factors  were  assigned  to  each  ring  of  injector 
holes  corresponding  to  the  relative  amount  of  energy  to  be  released  at 
that  location.  These  factors  are  shown  graphically  in  Fig.  11.  This 
factor  was  multiplied  by  the  number  of  holes  at  each  location  and  the 
orifices  for  each  location  were  selected  so  that 

w.  H.  F. 

—  .  (i  =  1-7)  (C-1) 

S  "i  ZHi  fi 

ttl 

where  w.  =  calibration  flow  rate  of  i —  orifice, 

1  ’ 

H.  =  number  of  holes  in  i —  location,  and 

F.  =  factor  of  i—  location. 

1 

The  values  of  flow  rate  for  each  orifice  and  the  arrangement  of 
orifice  combinations  for  each  profile  are  given  in  Tables  2  and  3  of 


Appendix  E. 
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APPENDIX  D 

METHOD  OF  EVALUATING  DATA 

The  output  signal  of  the  Photocon-Dynagage  system  was  simultaneously 
viewed  by  the  operator  on  an  oscilloscope  and  recorded  at  60  inches  per 
second  on  an  Ampex  tape  recorder.  After  the  run;  while  the  tape  was  re¬ 
played  at  1  7/8  inches  per  second,  the  signal  was  fed  into  the  Hathaway 
oscillographic  recorder.  %  comparison  of  the  incident  signal  with  a 
calibration  signal  of  known  frequency  and  amplitude,  the  frequency  and 
amplitude  of  the  pressure  oscillations  were  determined. 

Calculation  of  the  flow  rates  of  each  of  the  propellants  was  ac¬ 
complished  by  employing  the  following  equations: 

Wgtd  =0.668  (D-1) 

where  w  ,  ,  =  flow  rate, 

K  =  orifice  flow  constant, 

A  =  area  of  the  orifice, 

Y  =  orifice  expansion  factor, 

^  =  specific  weight, 

A  P  =  metering  orifice  pressure  drop 

and  suid  ),  are  based  on  an  upstream  temperature  of  520  R  and  the 

assumption  that  the  gas  follows  the  perfect  gas  law.  The  following  cor¬ 
rection  factor  was  used  to  correct  the  flow  rate  determined  by  equation  D-1. 
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(D-2) 


where  0,^  =  correction  factor, 

=  actual  upstream  temperature,  and 
C  =  compressibility  factor. 
Therefore, 


w 


^std 


(D-3) 


where  w  =  actual  flow  rate. 

The  equivalence  ratio  was  determined  by  dividing  the  actual  fuel- 
oxidizer  ratio  by  the  stoichiometric  fuel-oxidizer  ratio  for  the  propel¬ 
lant  combination. 

The  mean  chamber  pressure  was  obtained  directly  from  an  Esterline- 
Angus  Bourdon  tube  type  pressure  recorder.  Addition  of  atmospheric  pres¬ 


sure  converted  this  reading  from  gage  pressure  to  absolute  pressure. 


"Distribution  of  this  report  has  been  made  in  accordance  vith  the  Joint 
Army-Navy-Air  Force  Liquid  Propellant  Mailing  List  of  December  I961." 


